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FOREWORD 


“Thermal  Conduction  and  Ion  Temperatures  In.  the  Ionosphere" 
gives  the  results  of  a  study  of  the  ion  energy  balance  in  the  ionosphere. 
This  paper  will  be  published  during  1966  in  Barth  and  Planetary  Science 
latters. 


AVANT-PROPOS 


"Thermal  Conduction  and  Ion  Temperatures  in  the  Ionosphere" 
prdoenta  les  rdsulcats  d'une  dtude  du  bilan  dnergdtique  des  ions  dans 
l' ionosphere.  Ce  travail  sera  publid  en  1966  dans  Barth  and  Planetary 
Science  Letters. 


VOOPWOORD 


"Thermal  Conduction  and  Ion  Temperatures  in  the  Ionosphere" 
geeft  de  resultaten  van  een  ondersoek  betreffende  de  energiebalans  van 
de  lonen  in  de  ionosfeer.  Dit  artikel  cal  in  1966  in  Barth  and 
Planetary  Science  Letters  verschijnen. 


VORWORT 


"Thermal  Conduction  and  Ion  Temperatures  in  the  Ionosphere" 
glbt  die  Resultate  elner  Forschung  betreffs  der  Ioncnenergiebilans 
in  der  Ionosphere.  Diese  Arbeit  vird  im  laufenden  1966  in  Barth  and 
Planetary  Science  Letters  herausgegeben  verden. 


THERMAL  CONDUCTION  AND  ION  TEMPERATURES  IN  THE  IONOSPHERE 


by 

Peter  M.  BANKS 


Aba trace 


An  analysis  is  made  of  the  ion  energy  balance  equation  for  the 
ionosphere  including  the  effects  of  &  hM  and  H^u  ions,  thermal  conduc¬ 
tion,  and  the  cooling  by  atomic  hydrogen  and  helium.  It  is  shown  that 
thermal  conduction  strongly  affects  the  ion  temperature  profile  for  con¬ 
ditions  of  moderate  and  low  electron  densities  and  that  the  high  altitude 
ion  temperature  can  be  significantly  less  than  the  electron  temperature. 
Tn  addition,  1$  in  concentrations  as  low  as  10%  of  the  total  ion  density 
is  found  to  be  a  major  recipient  of  heat  from  the  ionospheric  electron 
°,a8  and  has  an  important  influence  upon  the  total  ion  energy  budget. 


RAsurad 


On  prdsente  une  analyse  du  biian  dnergdtique  des  ions  dans 

-f  f-  + 

1 ’ Ionosphere,  en  tenant  corapte  de  1 'influence  de  0  ,  Ea  et  H  ,  de  la 
conduction  therraique  et  do  re fro id is semen t  par  1 'hydrogen a  atomique  et 
par  I'hdlium.  On  montre  ainsi  quo  la  conduction  therm!  que  modifle  le 
profll  d;  la  temperature  ioniquo  pour  des  densitfis  Alectxoniques  moddrdes 
at  faibles.  A  haut2  altitude,  la  temperature  ionique  peut  8tre  nettement 
inf£rieure  &  la  temperature  dlectronique.  P' autre  part,  une  concentration 
H  atteignant  a  peine  10%  de  la  concentration  ionique  totalc  est  capable 
d: absorber  une  grande  quantity  de  chaleur  ct  influence  d&s  lore  la  biian 
dn&3‘gdcique  tota?  des  ions. 


2. 


Samenvafcting 


Er  wordt  een  analyse  gegeven  van  de  energlebalans  van  de  lonen 
to  da  ionosfeer.  Hen  houdt  hierbij  reksntng  met  de  Invloed  van  0  ,  He 
en  G  ,  met  de  therm!  *c he  geleiding  en  set  de  afkoeling  veroorzaakt  door 
ato«aire  water atof  en  helium.  Eveneens  wordt  de  sterke  invloed  aange- 
toond  van  de  thermlsche  geleiding  op  hut  ionentemperatuur-prof iol  voor 
gaaatigde  en  kleine  elektronendichtheden .  Op  prote  hoogte  kan  de  ionerj- 
temperatuur  juist  lager  liggen  dan  de  elektrcmentemperatuur .  Anderzijds 
is  de  concentratie  aan  H+,  die  nauwelljks  JOTS  van  de  totale  ionenccncen- 
tratle  bedra^gt,  in  staat  om  een  grote  hoeveelheid  warmte  op  te  slorpen  en 
aldua  de  totale  energiebalans  van  de  ionen  te  betnvloeden. 


Zusamoeofassung 


Die  Ionen  0  He  und  H  ,  die  WSraiaLeitung  und  die  Abktihlung 
durch  atomischen  Wasserstoff  und  Helium  werden  in  der  Ionenenergiebilanz 
berttcksichtigt.  So  sieht  man,  d&ss  die  Wttrmeleitung ,  ftfr  naesigen  und 
niedrigen  Elektronendicliten,  die  Ionentemperatur  verffndert  und  dass,  far 
grossen  HOhen,  die  Ionentemperatur  kieiner  als  dxe  Elektroncntemperatur 
aein  kann.  Anderseits,  flndet  man,  dass  eine  relative  Dichte  von  H+  von 
10  Prozent  eine  wichtige  QuantitKt  der  Wffrme  absorbiert  und  so  auf  die 
Bnergiebilanz  der  Ionen  einvlrkt. 


3. 


1.-  INTRODUCTION 


Through  numerous  experiment?  it  has  been  determined  that  the 
electron  temperature  above  150  km  is  consistently  higher  than  the  tempera¬ 
ture  of  the  neutral  atmosphere.  To  explain  this  state  of  thetmal  nonequi¬ 
librium,  a  number  of  theoretical  studies  have  been  made  of  the  electron 

heat  budget  equation  using  photolonization  as  the  principal  electron  gas 
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energy  source  L  *  '  '  '*  ■  .  The  most  recent  of  these  analyses’*  J  given 
results  in  general  agreement  with  current  measurements  but,  because  the 
electron  temperature  is  very  sensistlve  to  minor  changes  in  the  solar  EUV 
flux,  the  electron  density,  and  the  iiagnitude  of  the  assumed  mechanisms  of 
electron  energy  loss,  a  high  degree  of  correspondence  cannot  be  expected 
at  the  pruaent  time. 


The  primary  source  of  electron  energy  loss  above  300  km  altitude 

occurs  by  means  of  elastic  collisions  between  electrons  and  the  ambient 

ions.  The  ions,  in  turn,  transfer  this  excess  thermal  energy  to  the  neutral 

f3l 

gases  of  the  atmosphere.  Hanson*-  J  was  the  first  to  point  out  that,  because 
the  density  of  the  neutral  atmosphere  decreases much  more  rapidly  than  that  of 
the  electron  gas,  the  ion  temperature  should  have  a  pronounced  altitude  de¬ 
pendence  in  response  to  the  electron  heating,  increasing  gradually  from  the 
neutral  gas  temperature  near  300  km  to  the  electron  temperature  above  800  km. 


Experimental  support  for  the  transition  behavior  of  the  ion  temne- 

[6] 

rature  can  be  inferred  from  satellite  measurements  and  incoherent  rrdio 
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backscatter  data  L  It  is  necessary  to  point  out,  however,  that  although 

[3] 

the  results  of  Hanson1"  indicate  that  the  electron  and  ion  temperatures 
should  be  approximately  equal  at  high  altitudes,  in  fact,  in  the  available 
experimental  data*-  •  •  J  there  appear  to  be  significant  deviations  of  the 
ion  temperature  towards  values  which  are  consistently  lower  than  the  elec¬ 
tron  temperature,  even  in  the  atmospheric  regions  where  the  ion-neutral 
energy  losses  are  not  effective  in  cooling  the  ion  gas.  Further,  in  the 
regions  400-600  km  it  is  difficult  to  match  the  calculated  values  of  ion 


4.- 


t-esperatura  with  those  actually  measured  under  conditions  of  low  electron 
density. 


i'o  reconcile  the  differences  between  theory  and  experiment ,  a 
reanalysis  has  been  tade  of  the  ion  energy  balance  equations  to  Include 
the  previously  neglected  effects  of  ion  thermal  conduction,  the  presence 
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of  He  and  H  ions,  and  the  cooling  by  atomic  hydrogen  and  helium.  The 
results i  presented  in  this  paper,  indicate  that  Ion  heat  conduction  can 
play  an  Important  part  In  determining  the  profiles  of  Ion  temperature, 
leading  to  values  which  are  considerably  below  the  electron  temperature. 
It  Is  also  shown  that  H+  Ions  are  Important  In  the  over-all  ion  energy 
balance,  even  in  concentrations  less  than  10%  of  the  total  ion  density. 

II , -  ION  ENERGY  BALANCE 


The  time  dependent  energy  balance  equation  for  the  ion  gases  of 
the  ionosphere  is 


(1) 


where  U  -  3/2  n  k  T  is  the  ion  thermal  energy  per  unit  volume,  n  ■  n(0  )+ 

+  1  +  1  1  1 

n(He  )  +  n(H  )  ic  the  total  ion  number  density,  k  is  Boltzmann's  constant, 

T.  is  the  ion  temperature  (for  a  discussion  of  ion  temperature  coupling, 

■*  r10l 

see*'  ),  P,  is  the  rate  of  ion  gas  energy  production,  L.  is  the  rate  of 
energy  loss  to  the  neutral  atmosphere,  and  is  the  ion  heat  flux  arising 
from  conduction  and  diffusion  effects. 


Under  the  influence  of  thermal  conduction  and  the  non-local  hea¬ 
ting  effects  of  photoelectrons,  it  is  found  that  the  electron  tempe¬ 

rature  above  30C  km  is  significantly  greater  that  the  temperature  of  the 
neutral  atmosphere,  thus  providing  a  heat  source  for  the  ion  gases.  By 
means  of  electron- ion  collisions,  the  electron  gas  of  temperature  T  and 
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density  ne  transfers  energy  to  the  ions  at  the  rates, 

P. (0+)  -  4.8  x  10"7  n  n(0+)  [T  -  X .  ]  T  ’3/2  (2a) 

1  e  e  i  e 

P  (He")-  1.9  x  10‘6  n  n(He+)[T  -  T,]  X  '3/2  (2b) 

1  e  e  1  e 

P,(H+)  -  7.7  x  10'6  n  n(H+)  [x  -  X A  T  "3/2  ev  cnf3  sec'1,  (2c) 

i  e  e  1  e 

which  are  in  the  ratios  1:4:16.  Xhus,  the  heating  rates  fov  He+  and  H+  are 

equal  to  the  0+  rate  when  n(He+)/n(0+)  “  0.25  and  n(H+)/n(0+)  -  0.063, 

«■&» 

respectively,  Implying  the  electron  energy  transfer  to  the  He  and  H 
ions  plays  an  important  part  in  the  overall  ion  energy  balance. 

Xhe  rates  of  ion  energy  loss  tc  nrutral  gases  by  3»eans  of  elastic 
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collisions  have  been  discussed  recently  L  ’ ^  “  •vtd  these  results  arc  adopted 
for  this  study.  Briefly,  for  ions  in  parent  gases  the  process  of  resonance 
charge  exchange  is  dominant,  while  for  ions  in  unlike  gases  the  polarisation 
interaction  alone  is  assumed  te  be  acting.  Por  H+  and  0+  in  atomic  oxygen 
and  hydrogen,  the  effecto  of  accidentally  resonant  charge  exchange  have  been 
included. 


The  heat  flux,  q. ,  can  be  evaluated  for  a  magnetic  field-free 

1  r  13 1 

sedium  as,  neglecting  thermal  diffusion, 

qt  »-  K1VT1  +  5/2  n^T^  (3) 

where  i6  the  ion  thermal  conductivity  and  is  the  ion  diffusion  velo¬ 
city  relatte  to  the  mass  velocity  of  the  combined  ion  gas  and  the  neutral 
atmosphere.  The  evaluation  of  the  two  tenss  in  (3)  shows  that  heat  con¬ 
duction  is  generally  dominant,  but  that  there  are  some  circumstances,  espe¬ 
cially  for  low  temperatures  and  large  ion  densities,  when  the  diffusion  con¬ 
tribution  resulting  from  a  downward  flux  of  ions  cannot  he  neglected.  How¬ 
ever,  because  the  calculation  of  requires  the  solution  of  the  coupled  ic® 
temperature  and  density  equations  for  conditions  of  chemical  and  diffusive 
equilibrium,  it  is  difficult  to  adopt  reliable  vc.lues.  Hence,  in  this 
analysis  we  take  ■  0  and  deal  only  with  the  heat  flux  arising  from  ion 
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tcaperature  gradients  in  an  equilibrium  distribution  of  ion  gases. 

The  calculation  of  K  for  a  single  ion  gas  follows  from  the  work 
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of  Chapouu.  J  aa, 

K1  •«  4.6  x  104  T^2  /  k*/2  ev  cm"1  sec"1®^1  (4) 

where  Is  the  Ion  mass  in  atomic  units.  Unlike  the  electron  thermal  con¬ 
ductivity,  it  is  not  nacesaory  to  include  a  thermoelectric  reduction  factor 
or,  for  altitudes  above  300  km,  to  consider  the  effects  of  ion-neutral  col¬ 
lisions^1^. 


For  the  ion  gas  of  the  ionosphere,  an  extension  of  (4)  must  be 
made  to  include  the  effects  of  differing  conductivities  for  each  ion  species. 
An  exact  expression  is  difficult  to  obtain  and  for  the  present  analysis  an 
approximate,  density  weighted,  conductivity  has  been  used  in  the  form 

K,  »  1.15  x  104  T,5^2  [n(0+)  +  2n(He+)  +  4n(H*)]  /  n  ev  cm  1sec  1qK  1 .  (5) 

i  1  e 

When  the  effect  of  the  earth's  magnetic  field  is  considered  it  is 
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focndL  that  there  is  heat  flow  only  along  the  lines  of  geomagnetic  force. 
To  include  the  influence  of  the  inclination  of  the  field  lines  upon  ver¬ 
tical  profiles  of  ion  temperature,  we  Introduce  the  magnetic  dip  angle,  I, 
and  rewrite  (1)  as 


P  -  Li  +  8ln 


dT, 


dz 


thirl 


where  z  is  a  vertical  coordinate. 


(6) 


The  solution  of  (6)  for  the  ion  temperature  profile  requires  a 
knowledge  of  the  electron  and  ion  densities  in  the  upper  ionosphere.  Follo¬ 
wing  recent  work^10*1^,  we  adopt  a  state  of  0*  diffusive  equilibrium  above 
340  km  while  taking  H  and  He  to  be  in  chemical  equilibrium  with  N2»  02» 
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and  0  to  an  altitude  of  500  km,  above  which  a  condition  of  ternary  ion 
diffusive  equilibrium  is  applied" . 

The  numerical  solution  of  (6)  requires  a  choice  of  ion  temperature 
boundary  conditions  which  relate  to  the  Influences  of  ion  energy  production 
outside  the  region  of  analysis.  For  the  lower  boundary  condition  It  is 
assumed  that  conduction  is  ineffective  at  300  km  and  the  ion  temperature 
is  taken  as  the  solution  to  (6)  with  K  ■  0.  At  the  upper  limit 
(z  -  1500  km)  it  is  convenient  to  take  dT(/dz  ■  0  and  naglect  the  heat 

•a 

conducted  downwards  from  higher  regions. 

III.-  RESULTS 


The  ion  energy  balance  hao  been  solved  numerically  for  steady- 

r  is  i 

state  conditions  (ai^/dt  *  0)  using  three  atmospheric  models  of  NicoletL 
to  represent  the  effects  of  changing  aeronomic  conditions.  Throughout  all 
calculations  the  electron  temperature  has  been  treated  as  a  parameter  which 
has  the  same  value  for  all  altitudes. 

Several  points  of  Interest  have  been  found.  First,  it  appears  that 
the  ion  thermal  conduction  can  play  an  important  part  in  determining  the  ehape 
of  the  ion  temperature  profile  at  high  altitudes,  yielding,  for  a  constant  T  , 
isothermal  values  which  are  significantly  less  than  the  electron  temperature. 
This  effect  is  shown  in  curve  1  of  Figure  1  where  a  1000*K  neutral  atmosfrtiere 
has  been  used  with  only  0  ions  present.  The  electron  temperature  was  ta¬ 
ken  as  2500QK,  I  "  90°,  and  the  electron  density  determined  from  the  boun- 
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dary  condition  n  ■  2  x  10  cm  at  340  km.  For  curve  1  the  cooling  effects 
e 

of  atomic  hydrogen  and  helium  have  not  been  included.  To  show  the  impor¬ 
tance  of  thermal  conduction,  curve  2  (broken  line)  has  been  calculated 
with  -  0.  Since  this  choice  for  is  equivalent  to  the  conditions  found 
at  the  geomagnetic  equator,  it  is  seen  that  for  identical  aeronomic  para¬ 
meters  the  ion  and  electron  temperature  separation  would  nave  a  range  of 
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Fig.  I-  Ion  temperatures  with  an  electron  temperature  of  2500°K  and  neu- 

— -  tral  gas  temperature  of  1000°K,  for  a  magnetic  dip  angle  of  90° 

with  electron  denaity  of  ZxlO^cm"^  at  340  km.  Curve  1  inclu¬ 
des  (T  lens  and  thermal  conduction  but  neglects  cooling  by  He 
and  H*  In  curve  2  the  same  conditions  are  used  with  no  thermal 
conduction.  Curve  3  includes  0+  ions,  thermal  conduction,  He 
and  H.  Curve  4  represents  the  final  result  with  0+,  H+,  He,  H 
and  thermal  conduction. 
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350°1I  between  the  geomagnetic  equator  and  poles. 

A  second  effect  occurs  when  atomic  hydrogen  and  helium  are  ir<- 

+ 

eluded  in  the  0  energy  lose.  Curve  3  of  Figure  1  shows  that,  for  the  modal 
chosen  here,  the  high  altitude  cooling,  along  with  thermal  Conner  iron, 
reduces  the  isothermal  value  of  the  ion  temperature  to  1550*K.  The  intro¬ 
duction  of  h"*"  as  part  of  the  ionic  composition  reverses  this  situation > 
however,  by  providing  an  important  source  of  ion  heating  at  high  altitudes. 
As  shown  oy  curve  4  of  Figure  1,  this  results  in  an  increase  of  430°K  in 
the  isothermal  ion  temperature-  In  fact,  the  profile  obtained  for  the 
ion  temperature .omitting  the  effects  of  H+  is  very  similar  to  the  pro¬ 
file  found  when  H+  ,  He,  and  H  are  included.  The  basis  for  this  near 
equalicy  lies  in  the  counterbalancing  of  the  superior  heating  rate  of  H* 
by  the  large  H  charge  exchange  energy  loss  rate  to  atomic  hydrogen  and 
an  improved  ion  thermal  conductivity.  The  important  point  to  note  is  that 
for  the  adopted  model,  the  Ion  temperature  is  calculated  to  be  about  5C0°K 
less  than  the  electron  temperature.  Tba  effects  of  He  ,  which  are  net  shown 
here,  do  not  appear  to  be  large  and  can  be  neglected  in  a  first  analysis. 

The  value  of  the  electron  density  is  important  in  determining  the 
ion.  temperature  profile.  At  high  altitudes,  where  ion  cooling  by  the 
neutral  gases  is  ineffective,  conduction  must  carry  away  the  heat  obtained 
from  the  electron  gas.  The  changes  brought  about  by  varying  the  electron 
density  boundary  condition  at  340  km  are  shown  in  Figure  2  for  ng«  lxlO3, 
2xl05,  5xl05,  and  IxlO6  cm"3  with  T  -  1000°K,  T  -  2500°K,  and  both  0+  and 
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H  ions  being  included.  It  is  seen  that  ion  heat  conduction  is  generally 
important  during  conditions  of  low  electron  density  or,  more  exactly,  when 
the  high  altitude  energy  production  rate  is  small.  For  enchanced  electron 
densities  at  high  altitudes,  the  ion  conductivity  is  not  large  enough  to 
transport  downwards  the  energy  produced  in  the  ion  gac  following  moderate 
separations  between  the  ion  and  electron  temperatures. 
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B€f ectr,  of  changing  electron  density  boundary  conditions  upon  ion 
temperatures  with  the  electron  density  varying  from  10^  to  10^  cm 
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Changes  in  the  neutral  atmosphere  affect  the  ion  temperature  pro- 

Liles  by  shifting  the  region  of  ion-neutral  temperature  decoupling  with 

respect  to  the  electron  density  profile.  This  effect  is  shown  in  Figure  3 

where  modal  atmospheres  corresponding  to  neutral  thermospheric  temperatures 

of  .,00° ,  1000° ,  and  1394°K  have  been  used  with  T  -  2500®K  and  n  ■  2xl0^cm  ^ . 

e  e 

The  result,  which  is  due  to  a  combination  of  thermal  conduction  and  change 
in  the  ion  heating  rate,  is  that  for  given  values  of  electron  temperature 
and  density,  the  high  altitude  ion  temperature  decreases  slightly  with 
Increasing  neutral  gas  temperatures. 

As  a  final  point,  a  comparison  has  been  made  between  the  theore¬ 
tical  calculations  based  upon  (6)  with  0*  and  h"**  ions,  and  the  radar  ob¬ 
servations  of  Evans^'L  By  choosing  n  *  3.4  x  10"*  cm  ^  at  300  km, 

T  -  2600°K,  and  T  -  1000°K  to  match  the  observed  parameters  for  120C  EST, 
e 

September,  1963,  the  results  shown  in  Figure  4  have  been  obtained.  The 
squares  represent  the  experimental  data,  while  the  solid  curve  represents 
the  theoretical  result. 

To  evaluate  the  influence  of  diffusion  (see  discussion  following 

(3) ),  which  has  not  been  included  here,  the  divergence  of  the  second  term  In 

(4)  has  been  computed  using  the  final  results  of  the  ion  temperature  and 
electron  density  profiles.  For  the  diffusion  effect  to  have  an  importance 

equal  to  that  of  thermal  conduction,  it  is  found  that  the  ion  diffusion 

—  3-1 

velocity,  ,  would  have  to  be  about  7x10  cm  sec  which,  for  the  con¬ 
ditions  indicated  by  Evan's  data,  appears  to  be  too  large  to  represent  the 
true  physical  situation. 

As  a  test  for  the  existence  of  large  ion  diffusion  velocities  which 

could  alter  the  ion  temperature  profile,  it  is  noted  that  a  net  motion  of 
3  -1 

7  x  10  cm  sec  would  yield  a  doppler  shift  of  nearly  200  cycles/sec  at  a 
radar  frequency  of  430  Me/ sec.  Hence,  along  with  an  asymmetry  in  the 
characteristic  electron  temperature  peaks  there  should  occur  a  detectable 
shift  in  the  radar  spectra  half  power  points. 
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Fig.  4.-  Comparison  of  calculated  and  measured  ion  temperatures.  1’ne  so- 

- : -  lid  curve  was  calculated  using  a  1000°K  model  neutral  atmosphere 

and  the  experimental  parameters  of  Evam»L7J.  The  squares  indica¬ 
te  the  measured  values. 
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IV.-  CONCLUSIONS 


From  the  present  analysts  it  is  concluded  that  thermal  conduction 
within  the  ion  gases  of  the  ionosphere  can  be  important  in  maintaining  the 
high  altitude  ion  temperature  at  values  less  than  the  electron  temperature. 
Further,  an  ion  energy  budget  which  does  not  include  both  H  (  and  perhaps 
He+)  and  the  presence  of  atomic  hydrogen  and  helium  cannot  be  expected  to 
yield  temperature  prefiles  which  agree  with  expex imental  data.  The  effect 
of  ion  cooling  at  high  altitudes  by  means  of  downwards  diffusion,  although 
not  included  here,  may  also  be  important  under  particular  conditions  in 
further  Increasing  the  reparation  between  the  ion  and  electron  temperatures . 

Changes  in  the  parameters  of  the  charged  and  neutral  atmospheres 
affect  the  importance  of  thermal  conduction  in  the  ion  energy  balance. 
Increases  in  the  atmospheric  temperature  lead,  for  identical  electron  tem¬ 
peratures  and  charged  particle  densities,  to  a  decrease  in  the  high  alti¬ 
tude  ion  temperature.  Likewise,  increases  in  the  over-all  electron  densi¬ 
ty  can  be  expected  to  increase  the  ion  temperature  and  decrease  the  ion  ana 
electron  temperature  separation. 


Since  the  effectiveness  of  the  ion  heat  conduction  term  depends 
upon  the  dip  angle  of  the  earth's  magnetic  field,  the  separation  of  tem¬ 
peratures  at  high  altitudes  should  be  larger  for  higher  geomagnetic  lati- 
tudes.  Thus,  when  the  data  of  F&rleyu  J,  taken  at  the  geomagnetic  equator, 
is  compared  with  the  results  of  Evans^^,  and  Dcupnik  and  Nisbet^^,  there 
is  some  indication  that  there  are  increases  in  the  temperature  separations 
for  the  higher  latitude  sites. 
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